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Abstract. Asynchronous rotation in binary stars produces non-radial oscil- 
lations that are known to cause observable variability on orbital timescales. The 
horizontal perturbations of the surface velocity fields are referred to as "tidal 
flows" . In this paper we illustrate the manner in which tidal flows perturb 
the surface velocity field from that of uniform rotation, using a one-layer stel- 
lar model for the calculations. We justify the validity of this simplified model 
by the striking similarity between the photospheric absorption line-profiles it 
predicts and observational data of the binary system a Virginis. The velocity 
perturbations are used to compute the mechanical energy dissipation rates, E, 
due to the shearing flows for the case of a massive (50+28 M©) binary system 
having a moderately eccentric (e =0.3) orbit. The largest value of E around pe- 
riastron phases is found on the hemisphere facing the companion. However, at 
other orbital phases the maximum E may migrate towards the poles. Assuming 
that E plays a role in the mass-loss characteristics of massive binary systems, 
this suggests that peculiar binaries such as HD 5980 and 77 Carinae may have 
a highly non-sphcrically symmetric mass-loss distribution which, in addition, is 
time- variable. 



Asynchronous rotation 



A system is in synchronous rotation when the angular velocity of orbital motion 
f2 equals the angular velocity of the star's rotation loq. In eccentric orbits, the 
degree of synchronicity varies with orbital phase. We use periastron passage as 
the reference point for defining the synchronicity parameter (3 per = aJo/£l per = 
0.02Pv ro t(l — e) 3 ^ 2 Ri 1 (l+e) _1//2 , where e is the orbital eccentricity, the rotation 
velocity v ro t is given in km/s, the orbital period P is given in days, and the stellar 
equilibrium radius Ri is given in solar units. Non-radial oscillations are excited 
when at any orbital phase j3 7^1. The resulting surface velocity field consists of 
a pattern of perturbations superposed on the unperturbed stellar rotation field. 
We refer to the horizontal components of the velocity perturbations, AV V , as 
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"tidal flows" . A didactic manner of viewing the tidal flows is to imagine that 
due to the perturbations, some portions of the surface may be rotating slightly- 
faster than the average rotation rate while other portions may be rotating slightly 
slower. The size and distribution of these "portions" depend on the modes in 
which the star is oscillating. 

A huge body of research exists on the topic of tidal interactions, starting 
with Darwin (1879, 1880). References and very interesting discussions may be 
found in Eggleton et al. (1998) and Ferraz-Mello et al. (2008). However, the 
focus is generally on the long-term evolution of the orbital parameters. We in 
turn focus on the short-term effects, especially those that may influence the 
observational properties of stars on orbital timescales. 

Our method consists of computing the motion of a Lagrangian grid of sur- 
face elements distributed along a series of parallels (i.e., rings with different 
polar angles) covering the surface of a star of mass Mi as it is perturbed by 
its companion of mass M2. The main stellar body below the perturbed layer 
is assumed to be in uniform rotation. The equations of motion that are solved 
for the set of surface elements include the gravitational fields of Mi and M2, 
the Coriolis force, and the gas pressure. The motions of all surface elements are 
coupled through the viscous stresses included in the equations of motion. The 
kinetic energy of the tides may be dissipated through viscous shear, thus leading 
to an energy dissipation rate, E. The rate of energy dissipation per unit volume 
is given by the matrix product E = —P v : Vv, where P,, is the viscous part of 
the stress tensor and v is the velocity of a volume element with respect to the 
center of the star (for details see Moreno & Koenigsberger 1999; Moreno et al. 
2005; and Toledano et al. 2007). 

The benefits of our method are: 1) we make no a priori assumption regard- 
ing the mathematical formulation of the tidal flow structure since we derive the 
velocity field v from first principles; 2) the method is not limited to slow stellar 
rotation rates nor to small orbital eccentricities; and 3) it is computationally in- 
expensive and contains only two free parameters, the viscosity and the thickness 
of the surface layer. 

Surface Velocity Field and Energy Dissipation rates 

In Figure 1 we illustrate the characteristics of the perturbations of the azimuthal 
velocity field, A V v at periastron passage for a M1+M2 =50+28 M binary sys- 
tem with P=19.3 d, e=0.3 and Ri=21 R . Shown is the hemisphere containing 
the sub-binary longitude, ip =0. The corresponding values of AV^ on the equa- 
torial latitude are plotted in the right panel, illustrating the alternating positive 
and negative nature of this field. The largest perturbation in the direction of 
stellar rotation occurs slightly west of the sub-binary longitude. The largest 
perturbation in the opposite direction occurs at a longitude of ~60°. 

The shear energy dissipation rate E depends on the absolute value of the 
velocity perturbations, thus its azimuthal distribution has twice as many maxima 
and minima. A map of E corresponding to the velocity field displayed in Figure 
1 is in Figure 2 (left), with the corresponding plot for the equatorial latitude 
(right). Since E is a function of AV^,, energy dissipation rates are large wherever 
the absolute value of the azimuthal velocity is large. 
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Figure 1. Left: Distribution of the azimuthal velocity perturbations, AV Vl 
at periastron passage for the 50+28 M Q model discussed in the text. Color 
scale is such that white indicates maximum perturbation in the direction of 
rotation while black indicates maximum perturbation in the opposite direc- 
tion. The sub-binary longitude is the vertical line in the middle of the map; 
the rotation axis is tilted with respect to the plane of the paper by 55°. Right: 
Plot of AV V for the equatorial latitude, where positive values correspond to 
perturbations in the direction of rotation. 




Figure 2. Same as previous figure except that here we illustrate the distri- 
bution of the energy dissipation rate, E, and the color scale such that white 
indicates maximum E and black indicates minimum values. 
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It is interesting to follow the changes over the orbital cycle in the E distri- 
bution. A selection of orbital phases is illustrated in Figure 3 for the same binary 
system displayed in Figures 1 and 2. The maps are oriented so as to show the 
longitude range that would be visible to an observer at the given orbital phases. 
Shortly after periastron ((f) =0.10, bottom right), the pattern remains similar to 
that shown in Figure 2 (periastron), although the actual values of E are larger. 
As the star approaches apastron, the values rapidly decrease and the pattern of 
maximum E breaks down into smaller zones and migrates towards the poles. Af- 
ter apastron, the zones of large E begin to grow once again and migrate towards 
the equator. 




Figure 3. Orbital phase-dependence of the E distribution over the stellar 
surface for the same model illustrated in Figures 1 and 2 for phases with 
respect to periastron passage as follows: 0.10 (bottom right), 0.30 (top right) 
and 0.80 (top left) and 0.95 (bottom left). The maps are oriented so as to 
show the range in longitudes that would be visible to an external observer 
and color coding is the same as in Figure 2. 
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Photospheric line-profile variability 

A natural question that arises concerns the validity of the one-layer approxima- 
tion upon which our model is currently based. One way to address this issue 
is by comparing the model results with observations of photospheric absorption 
line-profile variability. 

Figure 4 illustrates the Si III 4552 A absorption line observed at four dif- 
ferent orbital phases in the B-type binary system a Virginis (Spica; P=4.01 
days, e=0.1). The observations were performed at the Canada France Hawaii 
(CFHT) 3.6m telescope with the ESPaDOns spectropolarimeter at a nominal 
spectral resolution of R=68,000 on 2008 March 20-28 (Harrington et al. 2009, in 
Preparation). The variability consists primarily of the presence of "bumps" on 
the line profile that systematically change their location over time, as was first 
studied in detail by Smith (1985a, 1985b) is clearly seen. Noteworthy also is the 
"boxy" shape of the absorption at some orbital phases. The dotted lines super- 
posed on the observations correspond to the theoretical line profiles computed 
with our model for the corresponding orbital phases. The model reproduces both 
the general shape of the observational data as well as the number and location of 
"bumps". Since this model is computed directly from first principles with only 
the single-layer approximation as an assumption and with very few free param- 
eters, the similarity between the computed and observed line profiles is striking. 
The reason why the one-layer approximation works so well in this case is that 
the dominant effect producing the line profile variability arises from the exter- 
nal gravitational perturbation, rather than from the star's internal oscillation 
modes. 




Figure 4. Comparison of observations with theoretical (dots) line profiles 
for 4 orbital phases in the binary systems a Virginis (Spica). Note that the 
one-layer model yields line profiles with a general shape that is similar to 
that of the observations, and that the number and location of the sub-peaks 
also coincide. The absorption line from the lower-mass companion is clearly 
visible in the left panel, but in the right panel it blends with the primary's 
absorption. Its location is indicated by the label "m 2 " and accounts for the 
larger discrepancy between the observed and theoretical profiles in the right 
panel. 
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Tidal flows and surface activity 

Although synchronous rotation is generally assumed for the analysis of a large 
variety of binary systems, the uncertainties in the stellar radius, rotational ve- 
locity and orbital inclination generally allow for significant departures from the 
assumed corrotating state. Furthermore, no star in an eccentric orbit is in syn- 
chronous rotation. Hence, the presence of tidal flows on the surface of binary 
stars is likely to be much more prevalent than generally assumed. 

Tidal flows produce a surface velocity field having gradients in both the 
azimuthal and polar directions. In addition, due to the radial dependence of 
tidal forces, the velocity amplitudes on the surface are larger than those in 
internal regions, hence also creating a rotation velocity gradient in the radial 
direction. All these gradients are strongly time-dependent. We suggest that the 
strong and rapid variations in these gradients may lead to a much larger degree 
of turbulence on the stellar surface than is present in synchronized binaries or 
single stars. It is also tempting to speculate that these effects could potentiate 
the appearance of localized magnetic fields (Ulmschneider & Musielak 1998). 

Tidal flows lead to energy dissipation due to the associated viscous shear. 
The additional energy in surface layers could conceivably alter the radial tem- 
perature gradient as well as produce a non-uniform effective temperature dis- 
tribution over the stellar surface. It is interesting to speculate on the possible 
effect that this additional energy source could produce on stellar winds. For 
example, since the tidal shear is not uniformly distributed over the stellar sur- 
face, the question arises as to whether E can contribute towards enhancing the 
mass-loss rate. Alternatively, localized magnetic fields or surface convective re- 
gions might tend to inhibit the radiatively-driven outflows at the base of the 
wind. Hence, in extreme cases where E is significant, our model would predict 
a highly non-spherically symmetric mass-loss rate. In addition, the asymmetry 
would be greatest near periastron passage, with a wind whose characteristics are 
dominated by the strong perturbations near the equatorial latitudes while near 
aperiastron the effects of the perturbation might be concentrated near the polar 
regions. 
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